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Abstract 

We present a preliminary result of the branching fraction for the B meson decay to the final 
state K + ir~ir + via an intermediate K*° resonance using the sample of approximately 23 million 
BB mesons produced at the T(4S) resonance with the BABAR detector at the PEP II e + e~ collider. 
The K*° was detected through the decay to the final state K + ir~. The result of this analysis is 
B(B + -» K*°ir + ) = (15.5 ± 3.4 ± 1.8) x 10~ 6 where the first error is statistical and the second is 
systematic. 
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1 Introduction 



Charmless hadronic B meson decays to pseudoscalar- vector final states have contributions from tree 
and penguin diagrams that are of similar magnitudes. It is expected that the interference between 
the tree and the penguin contributions may lead to direct CP violation. The measurement of 
the branching fractions of these charmless decays is an important first step in our understanding 
of the sizes of the various contributions. Recent measurements of the decays of B mesons to the 
pseudoscalar-pseudoscalar states involving pions and kaons hint at larger penguin contributions 
than expected ||]. 

In this paper we report a preliminary measurement of the branching fraction for the decay 
B+ _> K*°tt + where the K*° resonance is detected through its decay to the K + ir final state. 
Inclusion of charge conjugate states is assumed throughout the paper. 

2 The BABAR Detector and Dataset 

The data used in these analyses were collected with the BABAR detector at the PEP-II storage ring 
Q. The BABAR detector, described in detail elsewhere ||, consists of five active sub-detectors. 
Surrounding the beam-pipe is a silicon vertex tracker (SVT) to track particles of momentum less 
than ~120MeV/c and to provide precision measurements of the positions of charged particles of 
all momenta as they leave the interaction point. A beam-support tube surrounds the SVT. Out- 
side this is a 40-layer drift chamber (DCH), filled with an 80:20 helium-isobutane gas mixture to 
minimize multiple scattering, providing measurements of track momenta in a 1.5 T magnetic field. 
It also provides energy-loss measurements that contribute to charged particle identification. Sur- 
rounding the drift chamber is a novel detector of internally reflected Cherenkov radiation (DIRC) 
that provides charged hadron identification in the barrel region. The DIRC consists of quartz bars 
of refractive index ~1.42 in which Cherenkov light is produced by relativistic charged particles. 
The light is internally reflected and collected by an array of photomultiplier tubes, which enable 
Cherenkov rings to be reconstructed and associated with the charged tracks in the DCH, providing 
a measurement of particle velocities. Outside the DIRC is a CsI(Tl) electromagnetic calorimeter 
(EMC) which is used to detect photons and neutral hadrons, and to provide electron identification. 
The EMC is surrounded by a superconducting coil which provides the magnetic field for tracking. 
Outside the coil, the flux return is instrumented with resistive plate chambers interspersed with 
iron (IFR) for the identification of muons and long-lived neutral hadrons. 

The data sample used for the analyses contains approximately 23 million BB pairs, corre- 
sponding to 20.7 fb _1 taken near the T(4S) resonance. An additional 2.6 fb -1 of data were taken 
approximately 40MeV below the T(4S') resonance to validate the contribution to backgrounds 
resulting from e + e~ annihilation into light qq pairs. These data have all been processed with 
reconstruction software to determine the three-momenta and positions of charged tracks and the 
energies and positions of photons. 

3 Analysis Method 
3.1 Candidate Selection 

The B candidates are reconstructed from three charged tracks. Charged tracks are required to have 
a transverse momentum greater than 0.1 GeV/c, at least 12 hits in the DCH and to originate close 



8 



to the beam-spot. A good quality vertex of the three tracks is also required which is 96% efficient 
for signal. 

One of the tracks is required to pass a tight kaon selection based on Cherenkov angle information 
from the DIRC combined with energy-loss information from the DCH. The track that passes the 
kaon selection is assigned the kaon mass. The other two tracks, which are candidate pions, are 
required to fail the kaon selection and are assigned the pion mass. The kaon selection passes 79% 
of signal kaons and 96% of signal pions pass the pion identification. 

The B candidates are required to satisfy kinematic constraints appropriate for B mesons pro- 
duced by the decay of an T(4S). We use two kinematic variables: the energy-substituted B mass 

ttjes = \J ( s /2 + Po • Pb) 2 /Eq — p B , where the subscripts and B refer to the e + e~ system and 
the B candidate, respectively; and AE = E* B — y/s/2, where E* B is the B candidate energy in the 
center-of-mass frame. For signal events, the former has a value close to the B meson mass and 
the latter should be close to zero. In the Monte Carlo simulation of B + — ► K*°ir~ the AE and 
rriES distributions are Gaussian and have have resolutions of 18MeV and 2.4MeV/c 2 respectively. 
The requirements on these variables, | AE |< 60MeV and | tubs — 5279.5 |< 7.5MeV/c 2 , form 
the signal region in the ities-AE plane. The bounds of the signal region are over three times the 
Monte Carlo predicted resolutions so that uncertanties in the AE and ities distributions do not 
contribute significant systematic errors. 

The K*° candidate is formed by adding the 4-momenta of the kaon and oppositely charge pion 
candidates. A window of 0.816 — 0.976 GeV/ c 2 is allowed for the mass of the K + ir~ pair, ±1.6 times 
the width of the K*° resonance. For B decays to final states formed from a pseudoscalar and a 
vector meson, the vector meson is longitudinally polarized. The angular distribution of K*° decay 
products is used to suppress background. The helicity angle, Oh, is defined as the angle of the K 
momentum in the K*° rest frame with respect to the momentum of the K*° in the lab frame. The 
distribution of Oh follows a cos 2 Oh functional form for signal and is essentially flat for background. 
A requirement of | cos Oh |> 0.2 is used in this analysis. 

To avoid bias, the on-resonance events in the signal region were not counted or studied until the 
analysis method and selection criteria were finalised. Many of the selection criteria were optimized 
for statistical significance to an assumed branching fraction, as the signal region was blinded. The 
criteria that were optimized are the mass window on the K* candidates, cos Oh, the criteria on the 
particle identification and the requirements on event shape discussed in the following section. 

3.2 Background Suppression and Characterization 

There are two contributions to the background. The dominant background is random combinatorial 
background but there is also a small contribution from specific B decay channels. 

Charmless hadronic modes suffer large amounts of background from random combinations of 
tracks, mostly from light quark and charm continuum production. Such backgrounds may be 
reduced by selection requirements on the event topologies computed in the T(45) rest frame. As 
the B mesons are almost at rest in this frame, they decay spherically, whereas continuum events 
are jet-like. We use cosOt, the cosine of the angle 0t between the thrust axis of the B meson 
decay and the thrust axis of the rest of the event (all tracks and neutral clusters in the event except 
those that form the B candidate). For continuum-related backgrounds, these two directions tend 
to be aligned because the fake reconstructed B candidate daughters lie in the same jets as the 
remaining particles in the event. By contrast, in B events, the decay products from one B meson 
are independent of those in the other, making the distribution of this angle isotropic. A requirement 
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on the maximum size of | cos $t\ strongly suppresses continuum background. This requirement was 
optimised to give the greatest significance to the branching fraction. | cos Ot |< 0.7 was optimum 
which is 67% efficient for signal and rejects 90% of continuum events. 

Other event shape variables also help to discriminate between bb and continuum events. We 
form a linear combination of 11 variables in a Fisher discriminant T Q. The coefficients for each 
variable are chosen to maximize the separation between training samples of signal and background 
events. The variables contained in T are the cosine of the angle between the B momentum and 
the beam axis, the cosine of the angle between the thrust axis of the B candidate and the beam 
axis, and the summed momentum of the rest of the event in nine cones coaxial with the thrust axis 
of the B candidate. The J- requirement is 80% efficient on signal and rejects 60% of continuum 
events after the cos 9t cut has been applied. 

Despite the power of such topological variables to reduce the combinatorial backgrounds, it 
is necessary to make a background subtraction in measuring the branching fraction. In order to 
do this, the background in the signal region is estimated from the number of on-resonance data 
events in a sideband region and extrapolated into the signal region. The sideband region is a 
larger area than the signal region at lower tties and is specified by 5.22 < thes < 5.26GeV/c 2 , 
\AE\ < 0.3 GeV. In order to extrapolate from the sideband region to the signal region the shape of 
the background in the AE-ities plane must be characterised. This is done using on-resonance data 
(avoiding the signal region), off-resonance data and Monte Carlo continuum events in the region 
5.20 < niES < 5.29 GeV/c 2 , \AE\ < —0.3 GeV. The shape of the distribution of the background as a 
function of toes is parameterized using the ARGUS function Q which has one shape parameter, £. 
The background AE distribution is quadratic and the integrals over the signal region and sideband 
depend only on the quadratic coefficient. The quadratic coefficient and £ values measured for on- 
and off-resonance data and continuum Monte Carlo (which are consistent) are averaged and used 
along with the AE and rngg requirements to calculate 1Z, the ratio of the number of candidates in 
the signal region to the number in the sideband region. For this analysis, 1Z = 0.0647 ± 0.00029. 

In order not to underestimate the background, the selection criteria optimisation is done using 
only half the on-resonance events (every alternative event) and the final number of events in the 
sideband is measured on the other half of the events. 

The B background is supressed by the particle identification, and the K*, AE, and rriES require- 
ments. The remaining B background is then estimated with Monte Carlo events and subtracted 
from the number of events in the signal region. Possible sources of l?-related backgrounds include 
events with low-multiplicity decays to charm and other charmless decays. Background channels 
were identified by running over generic bb Monte Carlo and a dedicated charmless B decay simula- 
tion with loosened analysis requirements. For each background, Monte Carlo events exclusively of 
that type were studied to predict feed-through. Current branching fractions || are used to calculate 
the number of events expected in the signal region. Three channels have potential contributions, 
B + — > K*°K, B + — > p°K + and B + — > p°Tr + though they are all small. The total amount of 
-B-related background expected in the signal region is only 1.0 ± 0.6 events, the error is from the 
errors on the branching fractions and is systematic. 

Non-resonant B + — > K +, k~tt + decays and decays to a K + tt~it + final states proceeding through 
higher K + ir~ resonances are expected to be negligable in the signal region (less than one event). 
Although backgrounds from channels with the same final state are small the possible interference 
of the amplitudes of these processes may affect the branching fraction in the area of the Dalitz plot 
under study. This effect is not considered in the the branching fraction calculation but could be as 
large as 30%. With more data this interference can be studied. 
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3.3 Branching Fraction Analysis 



The branching fraction is calculated using 

B _ N 1 -2KN 2 -N Bbg 
N BS xe 

where Nx is the number of candidates in the signal region for on-resonance data; N2 is the number 
of candidates in half the on-resonance data observed in the sideband region, so that 2IZN2 is the 
estimated number of background candidates in the signal region; Nsbg is the number of 5-related 
background expected in the signal region; N Bg - is the number of BB pairs produced and e is the 
signal efficiency multiplied by B(K*° — > K + n~) = 2/3. 

For the signal efficiency in Eq. |], we used simulated signal events, applied the same selection 
criteria as used for data, and corrected for discrepancies between Monte Carlo and data. The 
efficiencies due to tracking, particle identification, vertex quality and T selection criteria are deter- 
mined from independent control samples derived from the data. 

Any discrepancy between the measured track-finding efficiency and that obtained in our simu- 
lation is taken into account with a momentum, azimuthal angle, and track multiplicity dependent 
weight per track. The ir^ and K particle identification efficiency was measured in data using a 
control sample of D + —* D°ir + ,D° — > K + ir~ . The efficiency was measured as a function of particle 
momentum and polar and azimuthal angle. 

The efficiency of the T requirement was measured using the control sample B + — > D°ir + , D° — > 
K + n~ which has the same final state and similar kinematics to the signal mode. The vertex quality 
requirement efficiency was not measured on this channel as the D° may fly before decaying. Instead 
the control sample B + — > J/ipK + , J ftp — > was used. 



4 Treatment of Experimental Uncertainties 

The systematic errors associated with branching fraction measurement arise from uncertainties in 
the background subtraction, in the overall signal efficiency and in the number of BB pairs produced 
as this is only known to 1.6%. 

The estimate of the combinatoric background is given as the product of the number of events 
in the sideband region and the factor 1Z. The number of events observed in the sideband has only 
a statistical error. The factor 1Z has a systematic error given by the error in the fitted ARGUS 
shape parameter, £. The contribution from the error in the AE characterization is negligible. 
A systematic uncertainty on 1Z is estimated from the range of values of £ measured in different 
data and different AE regions. The error on 1Z contributes a 2.5% systematic uncertainty to the 
branching fraction. The error on N^bg the number of -B-related events expected in the signal region, 
is from the errors on the branching ratios used in the calculation. This contributes 1.8% to the 
systematic error. 

Some uncertainty arises from the limited statistics of the signal Monte Carlo used to estimate the 
efficiency. The fractional error from this source is 2.3%. The accuracy of the simulation is subject 
to systematic uncertainties leading to uncertainties in the efficiencies of selection requirements. At 
present the dominant uncertainties are in the the efficiencies of the particle identification (6.4%), 
vertex quality cut (6.0%), T cut (4.1%) and tracking (3.6%). These systematic errors come from 
the methods used to to calculate the efficiencies in data, explain above. 

The uncertainty on the particle identification data efficiency is taken as 5.0% per kaon and 2.0% 
per pion. The particle identification efficiency errors for the two pions are added coherently but the 
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Figure 1: AE and tues distributions for B + — > K ir + ,K 
ground distributions used in the background subtraction. 



K + ir . The lines are the back- 



pion and kaon identification efficiencies have little correlation and so are added in quadrature. The 
uncertainty in the determination of the tracking efficiency in data leads to an uncertainty in the 
tracking efficiency of 1.2% per track, added coherently for all charged tracks in the B candidate. 

The overall systematic uncertainty is the sum in quadrature of the contributions from all sources 
giving a systematic error on the branching fraction of 11%. 

5 Physics Results and Conclusion 

Our preliminary measurement result of the branching fraction is summarized in Table El Figure |l| 
shows the AE and ities distributions in the signal region with lines showing the background 
distributions used in the background subtraction. Figure |2| shows the background subtracted K*° 
candidate mass and cos##. The lines on the plots are the expected B + — ► K*°tt + distributions. 
The data is consistent with being polarised K*° mesons as the x 2 P er degree of freedom is less than 
unity for both distributions. 



Events in signal region 
Events in sideband 
Ratio, 1Z 

Combinatoric background 
i?-related background 
signal yield 
signal efficiency 


55 
148 

0.0647 ±0.0029 
19.1 ± 1.6 ±0.9 

1.0 ±0.6 
34.8 ±7.6 ±1.1 
(0.149 ±0.004 ±0.013) x 2/3 


Stat. Signif. (a) 
B (xlCT 6 ) 


6.0 

15.5 ± 3.4 ± 1.8 



Table 1: A summary of the branching fraction measurement analysis for the decay mode B + — > 
K*°tt + . The factor of 2/3 in the signal efficiency is the branching fraction of K*° — > K + n~ . 

In summary, we have measured the branching fraction for B meson decay to the final state 



12 



'.7 



0.8 0.9 1 1.1 

m Klt (GeV/c 2 ) 



-1 -0.5 0.5 1 

cos9 H 



Figure 2: The mass of the K candidate and cos Oh for on-resonance data in the signal region 
with background subtracted. The lines are the expexted B + — > K*®it + distributions. 



K + tt 7r + via an intermediate K resonance with a statistical significance of 6. Oct. Our preliminary 
result is B(B + -► K*°ir + ) = (15.5 ± 3.4 ± 1.8) x 10~ 6 . 
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